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We discuss the prospects of detecting the sources of ultra-high energy (UHE) cosmic ray (CR) 
nuclei via their emission of cosmogenic 7-rays in the GeV to TeV energy range. These 7-rays re- 
sult from electromagnetic cascades initiated by high energy photons, electrons and positrons that 
are emitted by CRs during their propagation in the cosmic radiation background and are inde- 
pendent of the simultaneous emission of 7-rays in the vicinity of the source. The corresponding 
production power by UHE CR nuclei (with mass number A and charge Z) is dominated by pion 
photo-production (oc A) and Bethe-Heitler pair production (oc Z"^). We show that the cosmogenic 
7-ray signal from a single steady UHE CR source is typically more robust with respect to variations 
of the source composition and injection spectrum than the accompanying signal of cosmogenic neu- 
trinos. We study the diffuse emission from the sum of extragalactic CR sources as well as the point 
source emission of the closest sources. 
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I. INTRODUCTION 



The origin and chemical composition of UHE CRs {E > 10^* eV) is a long-standing enigma in CR physics 1. 2]. The 
spectrum at these energies is expected to be dominated by extragalactic sources that have not yet been unambiguously 
identified. The average mass number of the composition can be inferred directly from atmospheric CR showers by 
measuring the elongation rate distribution and comparing this to simulations. Recent results of the Pierre Auger 
collaboration [3j indicate a transition of UHE CRs within the energy range 10"'^* eV to 4 x lO'^^ eV from a light 
spectrum (consistent with protons) towards a heavier composition [4]. In contrast, the HiRes collaboration 5J finds 
a mass composition compatible with that of a proton-dominated spectrum [3]. 

Indirect evidence for the CR composition can come from various features seen in the spectrum. The "ankle" at 
about 3 X 10^^ eV seems to be a natural candidate for the transition between galactic and extragalactic CRs fT'-^, 
but a lower energy crossover at the "second knee" at about 5 x 10^^ eV has also been advocated for proton-dominated 
spectra [lOl[TT]. Proton-dominance beyond the ankle is expected to be limited beyond the Greisen-Zatspin-Kuz'min 
(GZK) cutoff [T21[T3] due to resonant pion photo-production in the cosmic microwave background (CMB). Intriguingly, 
a suppression of the CR spectrum at the expected energy of about 5 x 10"'^^ eV has been observed at a statistically 
significant level [3 [M] and is consistent with a proton dominance at these energies. However, a similar feature could 
also originate from nuclear photo-disintegration of UHE CR nuclei in the cosmic radiation background (CRB), or 
from an in situ energy cut-off of the injection spectrum. 

In light of these yet inconclusive and even controversial experimental findings it is important to consider alternative 
cosmic messengers associated with the production and propagation of UHE CR nuclei. The sources of UHE CRs are 
also expected to emit high energy radiation in the form of neutrinos and 7-rays due to their interaction with ambient 
matter, radiation or magnetic fields inside or in the vicinity of the acceleration region. The strength and spectral 
energy distribution of this emission depends, however, not only on the uncertain chemical composition and injection 
spectrum of UHE CRs but also on the source environment. This dependence on source parameters is absent for 
the fluxes of cosmogenic neutrinos and 7-rays that are produced by CR interactions with the CRB and intergalactic 
magnetic fields (IGMFs); these contributions can be directly related to the observed spectrum of UHE CRs. 

The dominant interactions of UHE CR nuclei (atomic number A and charge Z) with the CRB are nuclear photo- 
disintegation [T5Vll7j . pion photo-production [18 and Bethe-Heitler (BH) pair production [TH]. Photo-disintegration 
in the CMB leads to a fast break-up of heavy nuclei at energies of about A x 3 x 10^^ eV with an interaction length 
of the order of a few Mpc. The direct contribution of photo-disintegration to cosmogenic neutrinos and 7-rays is 
negligible [5D1 [5T] . The dominant channel for cosmogenic neutrinos is the photo- production of charged pions in the 
CMB and their subsequent decay [18]. This process can be approximated by treating the nucleons of the nuclei 
as free protons and neutrons with energy E/A resulting in an production threshold with CMB photons of about 
A X 5 X 10^^ eV. In the decay chain of the charged pion the three emerging neutrinos carry away about 3/4 of the 
total pion's energy. 
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Photo-nucleon interactions produce roughly equal number of charged and neutral pious. Gamma-rays and 
produced via tt" and tt* decay, respectively, subsequently cascade on CRB photons via repeated e+e" pair production 
and inverse Compton scattering. The net result is a pile up of 7-rays at GeV-TeV energies, just below the threshold 
for further pair production on the diffuse optical background. Another contribution to this electromagnetic cascade 
comes from the BH production of e^e^ pairs. The corresponding energy loss length of CR nuclei is minimal at 
energies of A x 2 x 10^^ eV and decreases with charge as Z"^ /A. 

The energy loss rate h = dE/dt in electrons, positrons and 7-rays determines the bolometric flux of the cosmogenic 
GeV-TeV 7-ray signal and cosmogenic neutrinos. The energy loss via pion photo-production at CR energy E can be 
approximated by the energy loss rate of nucleons N as hA,-fTT{E) ~ AbN^^T^{E / A). Since pion photo-production has 
a relatively high energy threshold the contribution of cosmogenic neutrinos and 7-rays from this channel typically 
depends on the maximal energy of the UHE CR emission. This can lead to a significant model dependence of the 
prediction. Bethe-Heitler e+e~ pairs are already produced at much lower energy. The interaction can be treated 
as a continuous energy loss with rate 6a,bh(-E') = bp,Bii{E/A), where &p,BH is the BH energy loss of protons. 
The Z^-dependence of this process and the low threshold makes this channel an important contributor to the 7-ray 
cascade. 

We will study in the following the relative contributions of BH pair production and pion photo-production to the 
7-ray signal for CR models involving heavy nuclei. In section|ll]we start with a brief review of the propagation of UHE 
CRs and the calculation of the spectra. The development of electromagnetic cascades and the corresponding diffuse 
7-ray signal for CR models are discussed in section |ni[ We will focus here on two CR models, a proton-dominated 
spectrum with a low crossover and an iron-dominated CR model motivated by the chemical composition inferred by 
Auger. In section |IV| we discuss the prospects to detect the closest sources of UHE CR nuclei as point sources of 
cosmogenic 7-rays. This possibility depends on the presence of an IGMFs during the development of the cascade as 
outlined in section Iv] We finally conclude in section VI 



II. PROPAGATION OF COSMIC RAY NUCLEI 



The main reactions of UHE CR nuclei during their cosmic evolution are nuclear photo-disintegration |15H17j , pion 
photo-production [TS] and BH pair production [TH] on CRB photons. The angular-averaged differential rate of a 
transition between nuclei of type i and j with energy Ei and Ej is defined as 

1 

-f,^,iz,E,,Ej) = ^ I dcos0 I de (1-/3 cos 0K(z,e)^^(e'), (1) 
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where e) is the energy distribution of isotropic background photons at redshift z and e' = £7(1 — PcosO) the 

photon's energy in the rest frame of the nucleus with Lorcntz boost 7 ~ Ei/Anip. Besides the contribution of 
the CMB we use the cosmic infrared/optical background (CIB) from Ref. 22J. Due to the cosmic evolution of the 
CRB density the interaction rates scale with redshift. Whereas the CMB evolution follows an adiabatic expansion, 
nj{z, e) = (1 + z)^ nj{0, e/(l -I- z)), we assume that the CIB evolution follows the star formation rate as described in 
the appendices of Ref. [23]. From Eq. (nj) we define the integrated interaction rate Ti^j{Ei) = J dEjji^j{z, Ei, Ej) 
and the total interaction rate ri{Ei) — zlj ^i^jiEt). 

Photo-disintegration of nuclei with large mass number A is dominated by the giant dipole resonance (GDR) with 
main branches A {A — 1) + N and A {A — 2) + 2N where N indicates a proton or neutron [T5Ul7j . The 
GDR peak in the rest frame of the nucleus lies at about 20 MeV for one-nucleon emission, corresponding to i^QQi^ — 
Ax 2 X e^Jv ^ 10^° GeV in the cosmic frame with photon energies e = EmeV meV. The secondary nuclei with atomic 
number A — 1 and A — 2 inherit the boost of the initial nucleon and lie close to the next GDR at E^~^ and E^~^, 
respectively. Hence, the initial flux of nuclei emitted from CR sources rapidly cascades down to lighter nuclei. This 
leads to a suppression of the flux above _Ei.cs,a/^- 

The most general evolution of primary and secondary nuclei in the CRB includes all possible photo-disintegration 
transitions between nuclides (A, Z) competing with the decay of unstable nuclides. For simplicity, we follow the work 
of Puget, Stecker & Bredekamp (PSB) jl6j and consider only one stable isotope per mass number A in the decay 
chain of ^*^Fe. At energies below 10 MeV in the rest frame of the nucleus there exist typically a number of discrete 
excitation levels that can become significant for low mass nuclei. Above 30 MeV, where the photon wavelength 
becomes smaller than the size of the nucleus, the photon can interact via substructures of the nucleus. Out of these 
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the interaction with quasi-deuterons is typicahy most dominant and forms a plateau of the cross section up to the pion 
production threshold at ^ 145 MeV. We use the reaction code TALYS ^24] to evaluate the cross sections aA^B of the 
exclusive processes (7, N), (7, 2N), (7, a), (7, Na) and (7, 2a) (N stands for p or n) for nuclides of the PSB-chain with 
10 < A < 56. For the cross sections of light nuclei with mass numbers A ~ 2,3,4 and 9 we use the parametrization 
provided in Ref. \25^. 

Resonant photo-nuclear interactions with CMB photons set in at energies of 5 x 10^^ eV per nucleon and becomes 
hence more important for low mass fragments of the photo-disintegration process. We follow the approach outlined in 
Ref. |25j and approximate the total photo-nucleus interaction by the isospin averaged interaction rate of free nucleons 
as TjH,-yTriz, E) ~ ATj\[^^^T^(z,E/A) j25| . We also assume that the participating nucleon is removed from the nucleus 
and regard this as a contribution to one-nucleon losses. The nucleon-photon interaction rates can be determined 
using the Monte Carle package SOPHIA We refer to the Appendices of Refs. [231 HZj for further details of the 

calculation. 

Another important energy loss of UHE CR nuclei is BH pair production via scattering off the CRB photons. Since 
this is a coherent process of the nucleons the energy loss scales as where Z is the charge number of the nucleus [19] . 
Bethe-Heitler pair production with differential rate 7^1, bh can be treated as a continuous energy loss process with a 
rate 

bA,BH{z, E) = J dE'{E ~ E')jA,Bii{z, E, E') , (2) 

where E and E' are the energies of the nucleus before and after scattering, respectively. The energy loss of nuclei 
can be related to the loss of protons as 6a,bh(£^) = Z'^bp,Bii{E / A). The energy loss length E/b{E) via BH pairs is 
typically much smaller than the interaction length of resonant photo-disintegration of heavy nuclei, which is of the 
order of (4/A) Mpc. Hence the primary nuclei will be fully disintegrated within a few Mpc. The light secondaries with 
small charge per nucleus, like protons and helium, do not exhibit a Z^-enhancement. However, photo-disintegration 
preserves the total number of nucleons. In particular, the full disintegration of an initial source spectrum of nuclei 
Qa{E) into protons is equivalent to a proton emission rate Qp{E) ~ A^Qa{EA). We will see later that the energy 
loss via BH of actual models of UHE CR nuclei receives also some -enhancement via BH loss for nearby sources. 

For the calculation of the diffuse spectra of UHE CR nuclei we assume that the cosmic source distribution is spatially 
homogeneous and isotropic. The comoving number density Yi = ni/{l + z)'^ of a nuclei of type i is then governed by 
a set of Boltzmann equations of the form: 

= dE{HEY,) + dEihY,) ~T,Y,+Y,j di^j 7j^<>S- + . (3) 

j 

together with the Friedman-Lemaitre equations describing the cosmic expansion rate H[z) as a function of the redshift 
z. This is given by H'^(z) = Hq [ilm(l + z)^ + ^a], normalized to its value today of Hq ^ 70 kms^^ Mpc""'^, in the 
usual "concordance model" dominated by a cosmological constant with ^l\ ^ 0.7 and a (cold) matter component, 
rim ~ 0.3 [2]. The time-dependence of the redshift can be expressed via dz — —dt (1 + z)H. 

The term Ci in Eq. ([s]) corresponds to the emission rate of CRs of type i per comoving volume. For the CR injection 
spectrum of nuclei of mass number A we use a power-law approximation with exponential cut-offs above i?max and 
below E^in, 

Qa{E) a E-^ exp{-E/E^,^) exp{~E^in/E) . (4) 

To account for cosmic evolution of the spectral emission rate per comoving volume we introduce an energy-independent 
scaling of the form Ca{z, E) — 'H{z)Q a{E) where we use the approximation 

H(z) = Ho(l + ZYQ{Z - Z,nin)e(Zmax ^ z) , (5) 

with Zmin = and Zmax = 2 unless otherwise stated. Note that our ansatz for Ca{z, E) does not distinguish between 
the cosmic evolution of the CR source density and the evolution of the intrinsic emission rate. This distinction is not 
important for the calculation of the diffuse CR spectra, but plays a role in the prediction of neutrino and 7-ray point 
source fluxes associated with these CR sources. We will assume that the emission rate of CR sources is fixed and that 
their number density evolves with redshift. 

In the following we are going to consider two models of extragalactic CR sources, that have been considered 
previously in fitting the UHE CR data [TOl HH] ■ The first model consists of CR proton sources with a strong evolution 
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FIG. 1: Left panel: Two models of extragalactic CRs assuming a homogenous distribution of protons (red line) and iron (blue 
line) between Zmin = and Zmax = 2. For the proton sources we use an injection spectrum with 7 = 2.3, i5min = lO^* eV, 
10^° '^ eV and assume strong source evolution with n = 5. The extragalactic iron sources assume an injection spectrum 
with 7 = 2.3, i?min = 10^* eV, i?max = 26 X 10^"'^ eV no evolution n = 0. Right panel: The corresponding spectra of 
cosmogenic 7-rays (dashed lines) and neutrinos (dotted line) for the two models. The diffuse 7-ray spectrum of the proton 
model is marginally consistent with the diffuse extragalactic spectrum inferred by Fermi-LAT [44] and the diffuse upper limit 
on cosmogenic neutrinos from the 40-string configuration (IC40) of IceCube [!?]• The cosmogenic 7-ray and neutrino spectra 
of the iron model are two orders of magnitude below the proton model predictions. 



(n = 5) with a relatively low crossover below the ankle. For the injection spectrum we use the power index 7 = 2.3 
and assume exponential cutoffs at iSmin = 10^^ eV and i?max = 10^°'^ eV (see Eq. (|4|). The spectrum of protons after 
propagation through the CRB is shown as a red line in the left panel of Fig. [T] The second model assumes a pure 
injection of iron with the same spectral index 7 = 2.3 but no evolution of the sources {n — Q). We assume the same 
exponential cutoff at low energies as in the case of the proton model, ii'min = 10^^ eV, and a high energy cutoff at 
-Emax = 26 X 10^°-^ eV, motivated by the rigidity dependence of the maximal energy of CR accelerators, i?max oc Z. 
The total spectrum of primary iron and secondary nuclei produced via photo-disintegration is shown as the blue line 
in the left panel of Fig. [T] 

Both models reproduce the UHE CR data above the ankle reasonably well. The deficit below the ankle is assumed 
to be supplemented by a galactic contribution. Note that the crossover with the galactic component is higher for 
the all-iron model than for the all-proton model. The fit of the model spectra to the CR data sets the absolute 
normalization of the CR emission rate. This can be expressed as the required bolometric power density per CR 
source, which depends on the local density of source, "Hq- For both models we find a value of 

-,42 ( '^0 \ ^ ^-1 



L ^ I dEEQiE) c 10- ( J^Z^^ ) ergs- . (6) 



III. ELECTROMAGNETIC CASCADES FROM HEAVY NUCLEI 



The evolution of cosmogenic electrons, positrons and 7-rays is governed by a set of Boltzmann equations analogous to 
Eqs. Elect romagnetic interactions of photons and leptons with the CRB can happen on time-scales much shorter 
than their production rates [29]. The driving processes of the electromagnetic cascade in the cosmic background 
photons are inverse Compton scattering (ICS) with CMB photons, -I- 7bgr — >■ -I- 7, and pair production (PP) 
with CMB and CIB radiation, 7-l-7bgr — >■ e+-|-e~ [TnKSD]- In particular, the spectral energy distribution of multi-TeV 
7-rays depends on the CIB background at low redshift. For our calculation we use the estimate of Franceschini et 
al. "SS]. We have little direct knowledge of the cosmic radio background. A theoretical estimate has been made by 
Protheroe & Biermann |31| of the intensity down to kHz frequencies, based on the observed luminosity function and 
radio spectra of normal galaxies and radio galaxies although there are large uncertainties in the assumed evolution. 
The calculated values are about a factor of ~ 2 above the measurements and to ensure maximal energy transfer in 
the cascade we will adopt this estimate and assume the same redshift scaling as the CIB. However, the 7-ray cascade 
below TeV does not significantly depend on the exact value of this contribution. A summary of the CRB used in this 
calculation can be found in Fig. A. 6 of Ref [35] ■ 
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High energetic electrons and positrons may also lose energy via synchrotron radiation in the IGMF strength B, the 
strength of which is limited to be below ~ 10~^G [33l|34j. Recently, the absence of (resolvable) GeV emission from 
TeV 7-ray blazars has been used to infer a lower limit on the IGMF strength of the order of lO^^^G The 
value of the IGMF has only little effect on the bolometric 7-ray flux in the GeV- TeV energy range relevant for our 
discussion [32j , but a sufficiently low value is crucial for a discussion of point-source emission. In the calculation of the 
cascade spectrum we will assume a weak IGMF with strength of 10^^^ G unless otherwise stated. Further processes 
contributing to the electromagnetic cascade are double pair production 7 -|- 7bgr + e~ + + e'^ and triple pair 

production, -I- jhgr —^e^+e~^ + e~. These contribution have only a minor effect on the 7-ray flux at GeV- TeV 
and are neglected in the calculation for simplicity. 

The observed diffuse extragalactic 7-ray flux is thought to be a superposition of various sources. Besides the 
cosmogenic 7-rays of UHE CR nuclei [ID] there are other candidates of truly diffusive processes associated with large- 
scale structure formation |41) or models utilizing the decay and annihilation of dark matter |42| . Other contributors 
are unresolved extragalactic 7-ray sources like active galactic nuclei, starburst galaxies, or 7-ray bursts (see Ref. |43j 
for a recent review). A recent analysis of the diffuse extragalactic 7-ray background (EGRB) by Fermi-LAT [M] 
shows a 7-ray spectrum that is lower and softer than previous results of EGRET [JS] . It has been argued that the 
Fermi-LAT flux constraints all-proton models of UHE CRs extending down to energies of the "second knee" |46) . 
though the systematics of UHE CR measurements is not sufficient to entirely exclude this model at a statistically 
significant level |32j . 

The right panel of Fig. [l] shows the diffuse 7-ray spectra (dashed lines) from the all-proton model (red lines) and 
the all-iron with (blue lines). In both cases we normalize the GR spectra (solid lines) to the Auger data above the 
ankle. The contributions from both models differ by about two orders of magnitude which is in qualitative agreement 
with the previous study [ID]. Whereas the all-iron model has only a negligible contribution to the EGRB the proton 
model saturates the observed background at 10-100 GeV. (In fact, decreasing the lower cutoff Ej^in = 10^* eV of the 
all-proton model would lead to an excess of the Fermi-LAT measurement.) For comparison, we also show cosmogenic 
neutrino flux (summed over flavors) of these models as dotted lines. The relative contributions from the two models 
differ by more than two orders of magnitude similar to the case of 7-rays. Note, that the cosmogenic neutrino flux 
from the all-proton model saturates a recent upper limit on the diffuse extragalactic neutrino ffux from the 40-string 
sub-array (IG40) of IceCube [47] . 

The strong model dependence of the diffuse fluxes is mainly due to the evolution of the sources as we will see in 
the following. The cascaded diffuse 7-ray ffux peaks in the GeV- TeV region and has an almost universal shape here. 
Its normalization can be determined by the total energy loss rate into 7-rays, electrons and positrons during the 
propagation of UHE GR nuclei. We can define the comoving energy density at redshift z as 

iOcUz) = J dEE[Y^{z,E) + Y,-{z,E) + Y,^{z,E)] , (7) 

which follows the evolution equation 

Wcas + ffWca. = J2 [ ^^iz, E)Ya{z, E) . (8) 

The energy density (eV cm~'^) of the electromagnetic background observed today is hence given by 



The relative effect of cosmic evolution on the energy density of the cascade can be estimated in the following way. 
The UHE CR interactions with background photons are rapid compared to cosmic time-scales. The energy threshold 
of these processes scale with redshift z as i?th/(l + z) where E'th is the (effective) threshold today. We can hence 
approximate the evolution of the energy density as 

Wcas+i?^cas-rycasH(^) j dE E Q{E) , (10) 

Bth/(l + z) 

where Tycas denotes the energy fraction of the CR luminosity converted to the electromagnetic cascade. Assuming a 
power-law injection Q(E) (x E"'^ with sufficiently large cutoff -Emax ^ -E-th we see that cosmic evolution enhances the 
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diffuse 7-spectrum as 

^cas(X J -^(l + z)"+7-3. (11) 


For the proton spectrum shown in Fig. [l] this corresponds to a relative factor of ^ 30. 

An additional, yet smaller relative factor depend on the chemical composition. We start with the the energy loss via 
photo- nucleon interactions &yi_^^(£') ~ Abf^^jr-f{E /A). This case is particularly simple to estimate: photo-disintegration 
losses conserve the total number of nucleons and we have the approximate relation A^Ya(£^A) ~ const. For the 
injection of a primary nucleus with mass number Aq and power-law index 7 we expect the scaling Wcas <x ■^o~^ for ^ 
universal high energy cutoff per nucleon. 

Energy loss by BH pair production follows the scaling 6a.bh(£') ~ Z'^bp,Bii{E /A). In the absence of photo- 
disintegration and photo-nucleon interactions this would result in a simple scaling of the form Wcas oc ZqAq~"'. 
However, as we have already discussed in section [IT] at those energies where pair-production is the dominant energy 
loss also nuclear photo-disintegration via the giant dipole resonance becomes important, shifting the average mass 
number and charge to lower values. On resonance, photo-disintegration has a typical inverse interaction rate of 
(4/A) Mpc and hence the primary nuclei will be fully disintegrated within a few Mpc. In this case we can expect that 
the dominant contribution to BH loss comes from light secondary nuclei and we would obtain the scaling Wcas oc Aq~'^ 
analogous to the case of photo-nuclear losses. 

For nearby sources the competition between photo-nuclear processes and BH loss makes it difficult to predict the 
exact scaling of these quantities. For distant sources and in particular for the calculation of diffuse spectra we expect 
that the scaling is closer to cjcas oc Aq~'' . For the diffuse 7-spectra of the iron model shown in the right panel of Fig. [l] 



this corresponds to a relative factor ^ 0.3 compared to the proton model. Together with the relative factor ( 11 ) from 
cosmic evolution and the difference in the normalization of the models (Qfc,GoV — O.OQp^Gcv) this accounts for an 
overall factor of ^ 110 in good agreement with the numerical value. 

In summary, the contribution of UHE CRs to the EGRB depends strongly on the underlying CR model, in particular, 
the evolution of the sources. The CR spectrum at the highest energy is, however, dominated by local sources. The 
relative contribution from these point-sources does not depend on the evolution of the full population. As we will see 
in the following, the predicted 7-ray flux from these source is relatively robust against model variation of the emission 
spectrum. 



IV. POINT-SOURCE FLUX 



So far, we have only considered diffuse 7-ray fluxes from a spatially homogenous and isotropic distribution of CR 
sources. However, the discreteness of CR sources can lead to local 7-ray excesses, in particular for the closest sources. 
Under optimal circumstances, i.e. sufficiently weak IGMFs, these excesses may even contribute as TeV 7-ray point- 
sources (PSs) in 7-ray observatories. This possibility has been previously studied for proton sources in Refs. |48H51j 
and has been revived recently in the context of unusually bright though distant TeV 7-ray sources [5^H5^ . Here we 
extend the discussion to the case of UHE CR nuclei and study the effect of IGMFs on the observability of the PS flux 
in detail. 

In the absence of an IGMF the flux from a PS at redshift z^, with emission rate Qa[E) [GeV~^ s~^] as in Eq. Q 
is equivalent to an integrated diffuse flux A'kJ{E) from a homogenous distribution on a sphere at comoving distance 
dc[z) = /g dz' /H{z'). The corresponding emission rate density is hence 

It is easy to check, that in the absence of interactions of the primary particle with the photon background the point 
source flux is then given as J*{E) = AttJ{E) = Qa{{^ + ^*)£')/(47ro?^). In particular, this reproduces the familiar 
luminosity-distance relation F — L/(47rd|) with flux F — J dEEJ*, luminosity L — J dEEQo and luminosity 
distance c?l = (1 + z)dc. 

For illustration, we show in Fig. [2] the 7-ray (solid lines) and neutrino (dashed lines) PS fluxes from a source at 
redshift z^, = 0.01 {dc{z^,) ~ 40 Mpc) emitting iron {A = 56), silicon {A = 28), nitrogen {A = 14) or protons. The 
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FIG. 2: The cosmogenic 7-ray and neutrino spectra of a CR point source of protons (red), nitrogen (green), silicon (yellow) 
and iron (blue) at a redshift z — 0.01 (~ 40Mpc). We assume an injection spectrum of the form Qa{E) oc E ' exp{—E/Emax) 
with index 7 and cutoff i5max as indicated in the plots. For comparison, we chose the same normalization of the nuclei in each 
plot such that L{E > 10^'' eV) = lO"^ erg/s (see Eq.je]). The solid line shows the total flux of 7-rays and the dashed line the 
cosmogenic neutrino flux. The dotted lines show the contribution to 7-rays from BH loss alone omitting the contribution from 
photo-pion interactions. 



different plots show variations of the spectral index 7 of the injection spectrum and the exponential cutoff i?max- 
In each plot the overall normalization of the proton emission rate Qp is chosen such that the source luminosity is 
Lp{E > lO^^eV) = erg/s (see Eq. [e]) and we use the same normalization constant for the other nuclei. The 
top panels show the results of 7 = 2.0 (top left panel) and 7 — 2.6 (top right panel) for a cutoff proportional to the 
nucleon mass -Emax = A x lO^*''^ eV. This choice corresponds to a universal exponential cutoff for the energy per 
nucleon in each injection spectrum. Since the energy loss Ba^'ytv from photo-nucleon interactions depend on E/A the 
spectra of cosmogenic neutrinos (dashed lines) have an almost universal shape as expected. Due to the universality 
of the neutrino spectra the relative normalization of their flux is in this case also given by Eq. ^ and scales as A'^~'^. 

The contribution of BH pair production to the 7-ray spectrum is shown separately in the plots as dotted lines. As 
discussed earlier, this contribution does in general not follow the behavior of the cosmogenic neutrino fluxes. 

If the initial nucleus is not fully photo-disintegrated, the BH contribution will be closer to Z^A^^''. Qualitatively, 
the variation of this contribution with spectral index 7 and initial mass number A is smaller than for the case of 
photo-pion loss and stays within a factor ~ 5. 

The robustness of this contribution becomes even more apparent in the case of a fixed maximal cutoff -Emax = 
2Q20.5 gY -[-Qj. g^jj nuclei which is shown in the lower panels of Fig. [2] for the same spectral indices 7 = 2 (bottom 
left panel) and 7 = 2.6 (bottom right panel). Since pion photo-production has a relatively high threshold of about 
A X 5 X 10^^ eV their contribution becomes strongly suppressed as we go to heavier nuclei and hence lower maximal 
energy per nucleon. This is apparent from the drastic decrease of the cosmogenic neutrino flux (dotted lines). 
(Similarly, the GZK 7-ray flux [551 156| at the upper end of the spectrum varies strongly with the maximal cutoff and 
compositon.) The 7-ray cascade, however, receives contributions from BH pair production at a lower energy threshold 



8 



[Mpc] Ho 


[10"'^ Mpc"^] 


L [10*° erg/s] 


4 


47 


2 


8 


6 


17 


16 


0.7 


137 



TABLE I: The local source density Ho (Eq. 141 and the source luminosity above lO'^^ eV (Eq. [gJ assuming a comoving distance 
Dj, to the closest CR source. 



and this contribution is only mildly effected by the variation of the model parameters. Generally, the comparatively 
small variation of the BH contribution in the electromagnetic cascade makes the prediction of 7-rays from the sources 
of UHE CRs more robust than cosmogenic neutrinos. 

The fit of UHE CR models to the data fixes the average luminosity density of CR sources. For the prediction of 
the average luminosity per source and hence the cosmogenic 7-ray flux from the closest CR source we have to fix the 
local source density Hq [cm~'^] introduced in Eq. The local density can not be much smaller than 10~^ Mpc~^ 
as can be estimated from the absence of "repeaters" in CR data [571 \SS\ ■ Moreover, the distant to the closest source 
can not be much larger than, say, 100 Mpc since UHE CRs will unlikely survive over longer distances. 

A spatially homogenous distribution of CR sources with number density 'H{z) per comoving volume as in Eq. ^ 
is equivalent to a diffuse flux of 

J{E)^j-Jdz^n{z)r{z,E), (13) 

where V — 47r(i^(z)/3 is the volume of the co-moving sphere containing the sources at redshift smaller than z. In the 
the following we will assume that the closest CR source is located at a comoving distance = dc{zi,). If the local 
source distribution with density Hq is sufficiently smooth we can expect that this is the only source within a distance 
Di = dcizi) > given by, 

no^^Jdz'-^H + zr^fDl (14) 







We can then decompose the diffuse flux as 



J{E)c^ ^J*{z,,E) + ^ f dz^n{z)J*{z,E). (15) 



Att An d 



\z 



We assume that the integral can be approximated by the diffuse flux of a spatially homogenous emission rate density 
L{z,E) = 'H{z)Qa{E) with Zmin — z\ as discussed in section [h] 

The emission of the closest source at D^, will only contribute to the UHE CR data at the upper end of the spectrum 
since 7^8. Due to the very poor statistic at these energies the position can only be determind within large statistical 
uncertainties. Instead - as our working hypothesis - we will assume in the following that the source location is fixed 
at Di, = Di/2. This corresponds to the average distance of a s ource uniformly distributed within a radius Di and 



weighted by the flux factor £)~^. With this choice the sum (15) will closely resemble the CR spectrum from a fully 
homogenous emission rate density with Zmin = as we will see in the following. We will study in the following three 
different source locations, D^, = 4, 8 and 16 Mpc. The corresponding source luminosities and volumes are tabulated 
in Table U 

The left panels of Fig. |3] show the contribution of the closest UHE CR sources (red lines) to the overall diffuse flux 

of UHE CRs for the case of the all-proton (top panel) and the all-iron (bottom panel) model. We also show as thin 

black lines the remaining contribution of the homogenous CR distribution beyond Di. Note, that in contrast to the 

PS proton spectra the total PS flux form the iron sources show a strong cutoff preceded by a small bump [59 . This 

limits the distance to the nearest UHE CR iron source to a few 10 Mpc in this model. The right panels of Fig. |3] show 

the corresponding PS fluxes in cosmogenic gamma rays (blue lines) and neutrinos (green lines). The PS 7-ray flux 

3/2 

from these close CR sources has the typical E^ ' form extending up to several tens of TeV following from inverse 
Compton emission of a fully Comptonized electron spectrum (E'^T^). 
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FIG. 3: Left panels: The average contribution of the closest proton (top panel) or iron (bottom panel) source to the spectrum 
of CRs. We assume that the closest source at determines the average local source density Ho as in Eq. [5] For illustration 
we assume local densities of 10""^, 10~* and 10~^ Mpc""^ with =4, 8 and 16 Mpc, respectively. From the fit to the CR data 
we can determine then the average source luminosity L. Right panels: The point-source 7-ray (blue) and neutrino (green) 
spectra for the closest proton (top panel) or iron (bottom panel) source. We assume again three different distances =4, 
8 and 16 Mpc. Whereas the diffuse 7-ray spectra shown in Fig. [T] of the two CR models differ by two orders of magnitude, 
the PS spectra are similar in magnitude for equidistant source locations. The PS fluxes of cosmogenic neutrinos depending on 
photo-pion interactions of protons are one order of magnitude lower for the case of an all iron source compared to an all proton 
source. 



The left panel of Fig.|4]shows the integrated 7-ray flux of these nearby CR sources in comparison with the sensitivity 
of present imaging atmospheric Cherenkov telescopes (lACTs) H.E.S.S. [30], MAGIC [SI] and VERITAS [SI], as well 
as the estimated future sensitivity of the water Cherenkov telescope HAWC [S3] and the Cherenkov Telescope Array 
(CTA) [64] ■ The solid and dashed lines show the 7-ray flux from proton and iron sources, respectively, at various 
distances. The signal depends only weakly on the compostion of the source. More important is the increased luminosity 
of the source (oc in the scenario of an increase local source density (cx D~^). For our three CR scenarios shown 
m Table |l] only a source distribution with a small local source density close to "Ho ^ 10"^ Mpc'^ (D* = 16 Mpc) and 
hence a large associated CR luminosity of L ~ 10**^ erg/s per source is expected to be visible in the future CTA after 
50h of observation. 

Note, however, that our ansatz D^, = Di/2 has been chosen for a good reproduction of the spatially homogenous 
emission density shown as the black lines in Fig. [3j It does not account for stochastic effects of the nearby source 
distribution. In general, the 7-ray signal of the closest CR source at comoving distance D^, and with fixed luminosity 
L is proportional to L/Di, and hence the model lines shown in Fig. |4] are expected to shift accordingly. Also, the 
detection of these multi-TeV 7-ray point sources in lACTs requires that the signal remains "point-like" , i. e. within 
the point-spread function (PSF) of the telescope. In this case it is important to consider the effect of an IGMFs on 
the development of the cascade as we will do in the next section. 
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E [GeV] [GeV] 

FIG. 4; Left panel: The sensitivity of present and future 7-ray observatories to nearby sources of UHE CRs for the all-proton 
and all-iron model. The integrated flux F of the point-source is almost independent of the the CR model considered. For a 
dilute UHE CR source density with Ho — 10'"' Mpc~'^ the future CTA should be able to identify the cosmogenic 7-ray flux 
of UHE CR sources within 50 hours of operation. Right panel: The size of the 7-ray halo for various IGMF strengths and 
7-ray energies. The halo size can be well approximated by a fit ^haio ~ 0.04°i3fG/-E7,TcV. 



V. EFFECT OF THE INTERGALACTIC MAGNETIC FIELD 



The cosmogenic 7-ray cascade of nearby CR sources can only contribute to a GeV-TeV PS flux if the deflections 
of secondary in the cascade via an IGMF is sufflcicntly small. We can estimate the extend of the cascaded 7-ray 
emission by simple geometric arguments following 36J. Deflection of electrons and positrons will be small if the 
energy loss length Ae of via inverse Compton scattering (ICS) is much smaller than the Larmor radius given 
as Rl = E/eB ~ l.l(£^Tcv/SfG)Mpc. (Here and in the following we use the abbreviations E — i^TcvTeV, etc.) 
For center of mass energies much lower than the electron mass, corresponding to energies below PeV in the CMB 
frame, electrons and positrons interact quickly on kpc scales but with low inelasticity proportional to their energy, 
Ae — 0.4 Mpc/i?TeV- The typical size of deflections of electrons and positrons is hence 9 ^ Xe/Rh ~ 0.2°i3fG/£'ioxev 

Deflection of close to the source have a smaller effect on the size of the halo then deflections close to the observer. 
To first order, if the cascade experiences a deflection /S.9 at a distance r from the observer, we can approximate the 
corresponding angular displacement A6'' in the observer's frame via Ad'/A9 ~ (d — r)/d. We can account for this 
scaling in the cascade equation by introducing the corresponding scaling in the Larmor radius R'j ~ RLd/{d — r) or, 
equivalently, by a scaling of the diffusion matrix of the form V ~ {{d — r)/d)^'D (see Appendix |A[). Since the distance 
of the closest CR source is expected to be smaller than the energy loss length by BH pair production or by photo-pion 
production the electrons and positrons will be produced continuously between the source and the observer. Hence, 
the average geometric suppression factor of the deflections is ~ 2 and hence (-R^) — 2Rl. 

In the GeV-TeV energy region the size of the 7-ray halo is almost independent of the source composition. This 
is a result of the rapid energy loss of via inverse Compton scattering above a few TeV compared to the slow 
production rate via BH loss of CRs or via pair production of 7-rays. The leptons quickly lose energy via ICS with 
CMB photons at a rate bjcg = E/X^; their spectrum in quasi-equilibrium {dtY^ ~ 0) follows the differential equation 
dEibics^e) — Pppi^- Thus, the Comptonized electron spectrum for E ^ -Emax has the form Ye ~ E^^. The deflection 
of an electron of the Comptonized spectrum is approximately 6e ^ Xs/Rl/^ following from dt{deYe) — Yc/Rl and 
Ye ^ E~'^ and (i?^) — 2i?i. The typical photon energy from ICS of a background photon with energy e is given by 
E^ ~ e{EelmeY and hence the halo is expected to extend up to an angle of about 0.01°i?fGemev/-£'7,TeV ~ independent 
of CR composition and source distance. 

We can define the halo size more rigorously with the approach outlined in Ref. [65] . We find that the halo size can 
be well approximated by the first moment of the angular distribution as 

0halo= \/2lfV>f'- (16) 

The result of the diffusion-cascade equation (see Appendix |A]) is shown in the right panel of Fig. |4j A fit to the 
data gives a numerical value of ^haio — 0.04°_BfG/-E7,ToV- This is consistent with our previous estimate for the most 




FIG. 5: Two models for the 7-ray spectra of the blazar source 1ES0229+200 located at redshift z = 0.14. The green data 
points show the H.E.S.S. observation and the green lines the estimated upper flux limits from the non-observation by Fermi-LAT 
inferred by Ref. [39 . Left panel: A model for the 7-ray spectrum assuming 7-ray emission at a rate oc i5~^''''0(2OTeV — -E). 
The solid green line shows the spectrum of secondary 7-rays without deflections in the IGMF. The dotted green lines indicate 
the part of the cascaded 7-ray spectrum within 0.1° around the source for an IGMF with coherence length As = 1 Mpc and 
strength Bq = lO^^'' G, 10"^^ G and lO"'^* G, respectively. Right panel: As in the left panel but now showing the 7-ray 
contribution from electromagnetic cascades assuming that the blazar is a CR proton source. We assume that the primary 7-ray 
emission is negligible and the (beamed) luminosity in protons is Lp ~ 10*^ erg/s. In this case the 7-ray flux is already below 
the Fermi-LAT upper limits and no IGMF is required to explain the data. 



abundant meV photons in the CMB spectrum. The typical size of the point-spread function (PSF) of lACTs is of the 
order of 0psF — 0.1°. Hence, an IGMF with a strength less than 10^^** G will not significantly decrease the sensitivity 
of future lACTs to the multi-TeV cosniogenic 7-ray signal of nearby CR sources. 

On the other hand, the non-observation of cascaded 7-rays as a GeV-TeV PS flux could imply a lower limit on the 
IGMF strength [23 . The observation of this effect requires that the sub-TeV 7-ray emission of the source is relatively 
quite, such that the cascaded spectrum would dominate the primary flux. Recently, the absence of (resolvable) GeV 
emission from TeV 7-ray blazars has been used to infer a lower limits on the IGMF strength at the level 10~^^G 

As an example, we consider here the emission of the blazar source 1ES0229+200 located at redshift z = 0.14, which 
has been detected by its TeV 7-ray emission by H.E.S.S. 66J. The spectrum is shown in Fig. [5] as the blue data. We 
consider two models for the 7-ray observation. In the right panel of Fig. [5] (from Ref. [55] ) we show a model assuming 
7-ray emission from the source with a rate oc i?^^/^0(2OTeV — E) (thin gray line). The surviving primary 7-rays 
are shown as a dashed green line and secondary cascaded 7-rays by a solid line. The cascaded spectrum would clearly 
dominate the sub- TeV emission and is inconsistent with upper limits from Fermi LAT (from Ref. [39, ). However, the 
signal within the PSF is significantly reduced by the presence of an IGMF as indicated by the dotted lines. 

The right panel of Fig. [Slshows an alternative model for the 7-ray emission of the blazar assuming strong emission 
of CR protons as in Eq. Ml with 7 = 2.3 and iS'max = lO^"'^ eV. This model has been advocated in Refs. [5^54j . 
In this case the observed spectrum is assumed to be dominated by cosmogenic 7-rays emitted during propagation. 
However, it is apparent that this model does not necessarily require the presence of an IGMF to be compatible with 
the Fermi-LAT limit as already noted by Ref. [SI]. We also show in this case the reduction of the PS signal via 
the presence of an IGMF with coherence length As = 1 Mpc and a strength Bp = 10^^® G, 10^^^ G and 10^^^ G, 
respectively. 

Note, that the required luminosity of the source is high in this case, ^ 10'*''^ erg/s, which is at least three 
orders of magnitude larger than the average luminosity of UHE CR sources inferred from the fit to the CR spectrum 
assuming a local source density larger than 10~^ Mpc^"^. In this model the blazar 1ES0229-I-200 can hence not be 
a typical source of UHE CRs. The CR emission along the blazar jet with opening angle 8 and the increase of the 
effective luminosity as 2/(1 — cos 5) does not play a role in this consideration since the same effect will also decrease 
the effective local source density of the anisotropically emitting CR sources. However, this example illustrates the 
strong model dependence on lower limits on the IGMF strength inferred by this method. 
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VI. CONCLUSION 



We have discussed the production of cosmogenic 7-rays in models of UHE CR nuclei. These 7-rays are a result of 
electromagnetic cascades in the CRB initiated by photo-pion production and BH pair production of the CR nuclei. 
The signal has its strongest contribution in the GeV to TeV range and is independent of CR interactions in the source 
environment prior to emission. We have discussed in detail how the 7-ray flux depend on injection spectra and the 
chemical composition of the sources. In general, we find that the flux of cosmogenic 7-rays is less model dependent 
than other agents of the CR interactions like cosmogenic neutrinos. 

As an illustration, we have studied two CR models of the UHE CR spectrum: an all-proton model with strong 
cosmic evolution and low transition to galactic CRs and an all-iron model dominating the spectrum beyond the 
ankle. The diffuse 7-ray flux from these models differs by two orders of magnitude. Whereas the proton model 
saturates the diffuse extragalactic 7-ray spectrum inferred by Fermi-LAT the iron model is practically unobservable 
in the background. We have shown that this large difference in the energy density of the cascade relies on the strong 
contribution of distant sources assumed in the all-proton model. 

The closest sources of UHE CRs can be observed via their 7-ray point source flux if the IGMF is sufficiently weak 
{Bq < 10"^'' G). We have argued that the 7-ray signal is expected to show only small variations with respect to the 
CR emission model due to the strong contribution of BH pair production. The absolute 7-ray flux depends on the 
CR luminosity and position of the source which can be related to the fit to the CR data and by estimates of the 
(average) local source density. We have estimated that the closest CR source should be observable via its cosmogenic 
7-ray emission in the future Cherenkov Telescope Array if the local source density is small {Ho ^ 10^^ Mpc^'') and 
hence the average source luminosity sufficiently large {L ^ 10^^ erg/s). 

We have also briefly commented on the possibility that the TeV emission of distant blazars can be naturally 
explained as a cosmogenic 7-ray signal if the blazar is a strong CR proton source (L ^ 10^^ erg/s). As an example we 
have studied the GeV- TeV emission of the blazar 1ES0229+200. The absence of strong GeV emission of this source 
has been used to derive lower limits on the IGMF strength. In contrast, the cosmogenic 7-ray emission of this source 
is consistent with the observation without the presence of an IGMF. 

We have only considered in this study steady sources of CRs, i.e. sources which have constant emission during the 
time of observation. Pulsed sources of UHE CRs could have a stronger 7-ray emission during the time of activity 
that could exceed our estimates. The sources of CRs are also expected to emit TeV 7-rays by CR interactions in 
the source environment. This contribution may dominate the point-source flux making the observation of cosmogenic 
7-rays difficult. However, in the case of pulsed CR sources, deflections of the cascade via magnetic fields can lead to a 
time-delay of cosmogenic 7-rays with respect to the in situ 7-ray emission of the source and may help to disentangle 
the contributions. 
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The small magnetic deflection of electrons and positrons in weak IGMF with coherence length and strength Bq 
can be treated as a diffusion process. We follow Ref. [5S] and deflne the moments of the 7-ray halo as 
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One can show that these moments folios the evolution equation 

= dE{HEYi''^) - T^Yi^\E) +J2 J d^'7/3c.(i?', + 4. J dE'V{E', E)Y^''-'\E') , (A2) 

/3=e,7_E E 

with diffusion matrix 

T>(F' p. , l min(l,ABrics(g)) e^i?g 

The quantity (x) denotes the inelasticity of of ICS and hence {x)ricg — 1/Aics- The first moment y'y^ equals the 
PS flux J*y^ as a solution of the Boltzmann equations (j3|. 

We define discrete values fJ"-* ~ AEiYe^\Ei), Q^ i ~ AEiQ^{Ei), etc. The combined effect of transitions and 
deflections within the cascade during a sufhciently small time-step At can be described by the matrix equations 

\YAt + At)J^ ^\T^,iAt) T,,{At)) ^\Y,{t)) ^ \Q-) 

(Y,it + At)Y'^y(T,,{At)n,iAt)\ (Y,it)Y\ ^JO 0\ (Y.it)^' ^ q) , (A5) 
The full cascade solution is then given by 

The 2n matrizes and S^") follow the recursive relation 

n 

{2PAt) = J2 -^^'^ (2^-^ At) • (2^-1 AO , (A7) 

n 

S(")(2PAi) = S(")(2P~iAt) + ^yl(*)(2P-iAt) • 6("-*)(2P-iAt) , (A8) 



1=0 



where the non-zero initial conditions are A^°\At) — T{At), A\j^ — diag{0, AW ij) and B^^\At) = 1. The matrices 

and B^^^ are the f amil iar tr ansfe r matrices for electromagnetic cascades in the presence of a source term. Using 
the recursion relations ( A7 ) and ( A8 ) we can efficiently calculate the matrices ^'"^ and yB*^"'' via matrix-doubling [57] . 
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